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Abstract 

Structure formation creates high temperature and density regions in the Universe that aUow the 
conversion of matter into more stable states, with a corresponding emission of relativistic matter and 
radiation. An example of such a mechanism is the supernova event, that releases relativistic neutrinos 
corresponding to 99% of the binding energy of remnant neutron star. We take this phenomena as 
a starting point for an assumption that similar processes could occur in the dark sector, where 
structure formation would generate a late time conversion of cold dark matter into a relativistic form 
of dark matter. We performed a phenomenological study about the limits of this conversion, where 
we assumed a transition profile that is a generalized version of the neutrino production in supernovae 
events. With this assumption, we obtained an interesting modification for the constraint over the 
cold dark matter density. We show that when comparing with the standard ACDM cosmology, there 
is no preference for conversion, although the best fit is within la from the standard model best fit. 
The methodology and the results obtained qualify this conversion hypothesis, from the large scale 
structure point of view, as a viable and interesting model to be tested in the future with small scale 
data, and mitigate discrepancies between observations at this scale and the pure cold dark matter 
model. 

1 Introduction 

Cosmological observations accumulated over the last decades have been favoring a picture of the Universe 
with 25% of the energy density made up of a weakly interacting non-baryonic particle. A scenario where 
cold dark matter (CDM) dominates is consistent with a bottom-up structure formation, therefore is 
favored by observations. This comes from the fact that the size of the initial perturbations is limited by 
the free-streaming length of the dark matter (DM) particles, for it is not possible to have gravitational 
clustering on scales larger then this length. 

The ACDM model with cold dark matter and a cosmological constant still provides the most simple 
description for the observations that nicely fits the data, e.g. the anisotropies of the Cosmic Microwave 
Background (CMB) measured by the Wilkinson Microwave Anisotropy Probe- 7yr (WMAP) [Tj, the 
matter power spectrum tracked by Luminous Red Galaxies (LRG) and measured by Sloan Digital 
Sky Survey-DR8 (SDSS) [2] and the recent accelerated expansion detected by the standard candles, 
supernova type la (SNIa), and collected by Union-2 compilation (Union-2) [3j. But, in spite of the 
ACDM success, there are still reasons to be cautious. A cosmological constant suffers from a serious 
fine-tuning and a coincidence problem [H El E] and, besides the dark energy puzzle, we haven't yet 
directly detected a dark matter particle, if there is any. Furthermore, specially regarding the small scales. 
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many controversies remain unsettled, e.g. the core-cusp [71 [8] and the missing satehites problem [9l [TO]. 
Generally, simulations of structures produce dark matter halos with a very steep profile, not matching 
observations, and the number of substructures around galaxies predicted by simulations is superior than 
the observed one. Efforts have been made in an attempt to improve both simulations and observations, 
but a general disagreement remains between approaches and results from different groups. 

Naturally, many alternative scenarios have been proposed to deal with the small scales issues, where 
dark matter, or at least a fraction of dark matter is not cold. The free streaming length of the DM particle 
limits the size of structures, setting the mass of a warmer candidate around the keV scale in a scenario 
in which this candidate is dominant. A dark matter particle with this mass scale could be a sterile 
neutrino jTTj. Some recent proposals of warm and warm+cold dark matter cosmologies have shown 
that a standard warm dark matter (WDM) cosmology does not solve the core-cusp problem p2j and, 
while WDM would naturally lead to a smaller number of satellites [13], other scenarios are not excluded, 
like the decay of cold thermal relics |14l I15j . On the other hand, the effect of baryons potentially plays 
a significant role in the galaxy dynamics [16] and still presents a challenge for simulations. Accounting 
for these effects can lower the number of subhalos as pointed out in ref . ( [171 [18] ) , however it is not yet 
certain whether a CDM+baryonic processes can reconcile simulations and observations of Milky Way 
satellites. More recently, there are claims |19j that gas outflows from stellar activity can significantly 
attenuate the steepness of the core. A proposal [20] of long range interactions between cold dark matter 
particles also would be able to dynamically settle the issue. Still in the spirit of mixed scenarios, it is 
interesting to see how cosmology can constrain limits on hot dark matter particles, like neutrinos or 
axions |21j . 

In this work we do not assess the particle nature of dark matter. We suppose that a cosmological 
scenario allows for a certain amount of hot dark matter and consider that this relativistic component 
is the result of a general conversion process. Naturally, the particles are not collisionless, but we don't 
make any restrictions about the interactions that could lead to conversion. Our work is inspired on pro- 
cesses in the baryonic sector that convert non-relativistic into relativistic matter. Known astrophysical 
phenomena can cause such conversions, like accretion disks around black holes and supernovae events, 
the latter one being the most interesting for us due to the large amount of energy released as relativistic 
neutrinos. 

When the gravitational pressure of a star exceeds the degeneracy pressure of electrons, the star 
explodes (or collapses) releasing radiation and neutrinos. The energy budget of the relativistic species 
is mostly carried by the neutrinos. Typically, this energy corresponds to 99% of the binding energy of 
the remnant neutron star. We take the net effect of this phenomena as a starting point for a supposition 
that similar processes can occur in the dark sector, and perform a phenomenological study about the 
limits in a late time and environmental conversion of cold dark matter into hot dark matter. There are 
many different models with such a transition, e.g. in ref. (\22\ I23j). where dark matter decays at earlier 
times and in ref. ( |24[ |25] ) where a coupling neutrino-quintessence accelerates neutrinos to relativistic 
velocities in scales of order of Mpc in the recent Universe. 

This paper is organized as follows. In section [2] we describe the model building. We present a 
model for the integrated diffuse neutrino flux in the Universe. Afterwards, we generalize the model 
for the conversion of cold to hot dark matter introducing the two parameters that describe the model: 
the initial time of the conversion Oj and the conversion rate rj. In section [3] we derive the Boltzmann 
equations that governs the dark matter particle's distribution, obtaining an equation for the energy 
density of both components. In section Owe describe the statistical analysis using Markov Chain Monte 
Carlo (MCMC) with the up to date data on large scale structure. Finally, in section [5] and [6] we discuss 
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our results and prospects. 



2 Supernova production rate 

Most models of supernova neutrino flux are built from direct observation of core collapse supernovae. 
However this method only provides a good parametrization for small redshifts, and since we are inter- 
ested in the effects of such production rate in all cosmological history, we will use a model where the time- 
dependent supernova rate RsNiz) is proportional to the star formation rate Rsf{z) p6l [271 128 } [29 l [30] 

„ , . 4>{m)dm 

^sn{z) = -:t25M^————Rsf{z) , (1) 
Jo.10 rn(f){m)dm 

and the star formation rate is parametrized as 

3.42 

RsFiz) = ^0-3 ^3.8,'^ ^^ [Mo.yr-^Mpc-3] . (2) 

We normalize the rate by considering the number of stars that would most likely undergo a supernova 
event with a correspondent emission of neutrinos, from 8Mq to 5OM0. This can be obtained trhough 
the integration of the initial mass function 

^i'^joc I ^_2.33_i.82iogm . o.lM0<m<lM0 • ^'^> 

where we used the combined Salpeter ^31j and Gould [32j initial mass functions in the appropriate mass 
intervals [33]. 

As it is well known, relativistic neutrinos are the main energetic output from a supernova explosion, 
where ~ 99% of the progenitor gravitational binding energy is transferred to neutrinos. As a conse- 
quence, a diffuse supernova neutrino flux is present in our Universe. The energy density carried by these 
neutrinos can be calculated from the evolution equations: 

Pb = -3Hpb - KRsN , (4) 
= -AHp^ + KRsN , (5) 

where H = d/a is the Hubble parameter and Ei, ~ 3 x 10^^ erg is the total energy carried by the 
neutrino flux in a supernova explosion. The solutions to these equations are given by: 



Pb 



1 
1 



/a , _ 
a^EyRsNdt 



(6) 



a 



p^ = — j a E^RsNdt . (7) 

In figure [T] we can see the profile of energy density carried by relativistic neutrinos produced by SN 
explosion as a function of redshift. In the same plot we show the energy density carried by relativistic 
species created in dark-matter decay, where the parameters were chosen to produce the same amount 
of energy density of relativistic species. It is possible to see from figure [1] that two very different 
mechanisms can produce similar patterns of relativistic species production. We propose a scenario 
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Figure 1: The supernova produced neutrino density fraction, Q = ph/pcr, as a function of z. 



where cold dark matter is converted into hot dark matter with a not initially populated distribution by 
a generic mechanism that is not necessarily the decay of the original cold dark matter particle content, 
but resembles the production mechanism of supernova neutrino. We chose the following ansatz for the 
cold dark matter density evolution 



Pea 



Poe 



-ri{a-ai) 



;a < Ui 
;a > Ui 



(8) 



where defines a starting time for the conversion and r], the rate that this conversion occurs. The 
evolution equations can then be easily integrated, and we obtain for the relativistic dark matter density 
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-r]{a-ai) 



;a < Qi 



(9) 



With this parametrization the conversion profile on SN neutrinos can be fairly reproduced by the the 
choice Oj = 0.3 and rj = 4.6 x 10~^, indicating that our parametrization presents a good versatility 
and a close connection to known astrophysical processes where a conversion from non-relativistic to 
relativistic matter occurs. 

When constraining our parameters, it is more convenient to work with variables that gives us directly 
the amount of dark matter in a relativistic form today. This can be calculated explicitly by: 



/ 



Ph 



Pc + Ph 



(10) 



a=l 



The correspondence between the parameter r] and /, for a given value of a^, is unique, and can be 
obtained numerically from the relations above. We present this correspondence in figure [2j 
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Figure 2: Values of parameter rj from / and 



3 Boltzmann equations 

Perturbations over the homogeneous background of the cosmic ingredients can be evolved as a series of 
coupled Boltzmann equations of perfect fluids, where each component is fully described as a hierarchized 
set of equations, truncated according to its thermal condition and coupled according to the physics 
considered. For a complete description of all the components, see the ref. ([M]). In this work, it is going 
to be detailed only the set of equations of cold and hot dark matter and their environmental conversion 
as a coupling. We are supposing the cold dark matter is a particle with mass rric that undergoes an 
environmental conversion at a rate dgdda to a hot dark matter particle with mass ruh with an empty 
distribution populated at a rate dgh/da. The Boltzmann equation for the cold dark matter is 

f=C[/J^ (11) 

In this scenario with dynamical dark matter, the total distribution function depends on the momentum 
and the scale factor, we also consider that the dependencies are detachable 

C[Up,a)] = -dUp)^^, (12) 

da 

where the function gc{o) represents an effective and already normalized conversion process. We consider 
that the phase space of the daughter particle is completely available. Integrating the distribution 
function, the final equations of density will 

. d .dgc{a) 

Pc + 3-pc = -aa— — pc, (13) 
a da 



and density perturbations are 

(5c + ikvc + 3$ = -aa^^4^5c. (14) 

da 

For non-relativistic particles is enough to extend the hierarchical description to the velocity term 

d ^ . dgJa) ^, 

Vc -\ — Vc + ik^ = —aa — - — Vc, (15) 
a da 

any higher order as pressure or stress is negligible and even undesirable, since there is no seen pressure 
effects in the clustering of dark matter. In a similar way, wc detail the equations of density and density 
perturbations of the hot dark matter, which is identical to the neutrino equations except for a non-null 
collision term 

^ = C[A]. (16) 

The distribution of relativistic particles in thermodynamical equilibrium is completely described by the 
temperature. We assume that this is the case and include the perturbations directly in the distribution 
function 

fh{x\pj,T) = fho{p,T) X [1 + ^{x\p,nj,T)]. (17) 
The collision term populates the hot dark matter from the cold dark matter distribution 

C[fhip,r)\ = -afcip) = afcip) . (18) 

The density of the relativistic particles will be given by 

d dgc (a-) 

Ph + 4:-Ph = aa— — Ph. (19) 
a da 

By integrating again the perturbed distribution function, but this time taking the expansion in Lcgendre 
polynomials, a hierarchical system of equations describing the perturbations in the hot dark matter 
temperature is obtained, where the monopole term is 

k-a%^^c-#V (20) 



and the dipole 

0h 



3 da 



(21) 



which can be used to obtain the whole set of equations by the recurrence formula 



A = [l^l-i -{1 + l)^+i] J > 2, (22) 

until some order of truncation Z^aa;) where 

''Tiiax J- h>n- ^max ''max -L V / 

For the desired density conversion given by equation 9, the ad-hoc conversion rate must be 

da a' 
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Table 1: Parameters of the model and their flat prior. The first six parameters in the first block belong 
to the so called "vanilla" cosmological model. The middle block contains extended standard parameters 
that are expected to have degeneracy with the conversion model parameters. The last block contains 
the parameters added by our model of dark matter conversion. 



Type 


Parameter 


Description 


Min 


Max 


\Aq Ti 1 1 1 a 




rvn Ti npiimfv 


0.005 


0.1 




n /?2 




n 04 


18 




e 


Ratio between the sound horizon and the angu- 


0.5 


10 






lar diameter distance at decoupling 








T 


Reionization optical depth 


0.01 


0.8 




Us 


Spectral index of primordial power spectrum at 
k = 0.05/iMpc-^ 


0.5 


1.5 






Amplitude of primordial power spectrum at A; = 
0.05/iMpc-^ 


2.7 


4.0 


Extended 


/. 


Fraction of neutrino density related to the dark 
matter 





0.1 




w 


Parameter of dark energy constant equation of 
state 


-1.5 


0.5 


Extra 




Scale factor when the conversion starts 


0.1 


1.00 




V 


Rate of conversion 





2 



which replaced in eqs. (|T3|) and p9j) . gives 



(25) 



Ph 



PcoW + V 



-1^ 



1 - e 



-ri{a-ai) 



(26) 



as planned beforehand in eqs. ([8]) and ([9]). Given that the conversion is supposed to be an environmental 
effect related to the galaxy halo, we include a scale dependence of the form 



da 
dghia, k) 



[l-h{k,kg)] X 



da 
dgh{a) 



(27) 



(28) 



da da 
where the term included is a smooth step function and kg is defined as the scale where the conversion 
reaches half of its maximum effect for each redshift. 



4 Statistical Analysis 

The statistical analysis was made by Markov Chain Monte Carlo (MCMC) using the data from the 
anisotropics of the cosmic microwave background measured by the WMAP, the matter power spectrum 
measured by the SDSS and the type la supernovae luminosity compiled by the Union-2. Best fits were 
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obtained by conventional analysis using Bayes posterior probability description. All priors used 
were flat distributions over the range tested. We run the CosmoMC |35j package of MCMC to test 
the theoretical predictions generated by a modified version of C AMB |36) , which included the particles 
described by the Boltzmann equations developed in section [3l Besides the standard parameters for 
cosmology, our model add the parameters listed in the last block of table [TJ The scale conversion kg 
was fixed to the scale of typical galactic halo size, kg = 0.2 Mpc-^ The convergence requirement 
followed the Gelman and Rubin R statistic, which states that the variance of chains means divided by 
the mean of chains variances must approach of one. In this work, we reached a convergence of at least 
R—l< 0.05, which could be regarded as a low convergence, but taken into account that the CosmoMC 
package is adjusted for statistical efficiency when fitting the standard cosmological model, we consider 
that for a first time fitting, the condition R — 1 < 0.05 for an alternative model is acceptable. 

Because there is an obvious degeneracy between neutrino density and any other kind of hot dark 
matter, we included the fraction of dark matter in massive neutrinos among the standard parameters 
of ACDM, distributing the density of neutrinos equally among three states and fixing the number of 
families to N^q = 3.04. It is also expected to exist a correlation between the effects of dark matter 
conversion and the parameter w of the dark energy equation of state [37], for this reason this parameter 
was included in our analysis. 

5 Results and Discussion 

By fitting the data with theoretical predictions, we calculate chains of the likelihood for all parameter 
values inside the ranges delimited. The complete dataset fitted was WMAP + SDSS + SNIa and the 
likelihood functions adopted for each observable were the ones suggested in their respective papers. We 
present in table [2] for each parameter p, the total mean value {{p)) of its distribution and its confidence 
level corresponding to the 95% central credible interval to represent the data. In tabled we also present 
the best fit value (j5) of each parameter and its confidence level corresponding to the 2a probable region. 
In this case, all parameters values are kept in the total best fit point. Limits not shown mean that the 
limit is compatible with the correspondent limit in its prior range. See ref. ([38]) for different statistical 
approaches and their relevance in cosmology. In order to compare and qualify the modified model, we 
also present the mean and best fit points, the credible regions and confidence levels of the standard 
cosmological model. We use AHCDM as the acronym for the model with dark matter conversion, while 
keeping ACDM for the standard model, both extended with massive neutrinos {rriu) and an arbitrary 
value for the dark energy equation of state parameter (w). When comparing the best fits, the minimal 
for the model with dark matter conversion depreciates the fit for the following amount 

= (XAHCDM - xXcDm) /^d.o.f. = 1.01/2. (29) 

For two extra degrees of freedom, this number means that the best fit for AHCDM is worse than the 
ACDM best fit, but it is withing la of confidence around the standard model best fit. The best fit point 
in the likelihood for the two extra parameters that characterize our model are 

= 0.6 , ai = 0.6 , (30) 

where the large standard deviation prevents better precision. These values for the global best fit 
correspond to a rate of converted hot dark matter over cold dark matter of around 50% in our Universe 
today. 
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Table 2: Mean total values (p) with their 95% central credible intervals and the best fits p with their 
confidence levels at 2(T, for the model with conversion (AHCDM) and for the standard model (ACDM). 
The fits were obtained with WMAP+SDSS+SN data. 



Type 


Parameter 


{p){95% CCI) 


p (2cj CL) 








ACDM 


Vanilla 




0222R-R^?f 

u.u^z,^Q 02X1 


0226^-^^?? 

u.uz,z,Uo.o214 


0223R-R^^§ 

U.Uz^^Og 0202 


0226^-^^;^^ 

u.Uz:^Uo 0209 


O /)2 


'-'•^'-'0.06 


'-'•-'-^"-'0.132 


'-'•'-"0.03 


u.i'±^0.122 


a 

u 


i.uoox 030 


1 033}-R^§ 

i.uoox 028 


i.uool 024 


i.uo'ii 024 


/ 


'-'•'-'■='0.06 


'-'•'-'"0.07 


'-'•'-'^0.04 


'-'•'-'^0.05 






'-'•^^0.96 


'-'•^<-'0.92 


-'-•'-"-'0.94 




3 llU^ 


q 1 q3.19 
"^•-■-"33.07 


"J- -'-'-'2.98 


"J- -"-^3.04 


Extended 


f 


0.03°-°^ 


O.O2O-O6 


O.O2O-IO 


O.O2O-O7 


W 


'-'•""1.06 


'-'•""-'0.98 


'-'•'-'"1.31 


'-'•"^1.13 


Derived 




68n«o 
"•""0.62 


66n«i 
"•""0.61 


70nI5 
"• ' "0.54 


68nli 
"•""0.58 


A P"p / r^vr 


1 o ol3.6 
^'^•"Jl3.0 


i2.6{i:i 


13 9I3.8 
iJ.Zx2.8 


12. 5x2:2 


0-8 


'-'•' J^O.63 


7fi0-85 
'-'•'O0.68 


'^0.58 


n 780.92 
'00.62 




11 14 


11 14 


11 16 




Ho 


69i 


' J-67 


7177 
' J-61 


7077 
'^65 




O.4O.8 


0.30-8 


O.O9I-4 


O.2I-I 


Extra 






N/A 


0.6 


N/A 




O.4O.8 


N/A 


0.61-0 


N/A 


Extra derived 




0.060-10 


N/A 


O.O7O-13 


N/A 


f 


O.4O.6 


N/A 


O.5O-7 


N/A 



In figure [3] we show the matter power spectrum from SDSS, the predictions for ACDM and the 
predictions for our model with dark matter conversion. In figure H] we show the correlation between the 
parameters Oj and r/. The confidence region for the two parameters combined was obtained while fixing 
the value of standard parameters to their total mean. In this case the set {oj, ry} excludes non- vanishing 
conversion {oj,??} — )• {1,0} at 2a. However, one should not take this as a preference for conversion. 
The other parameters are marginalized over in their total mean, which are significantly different from 
the standard model mean, therefore, the conversion is needed to restore the desirable evolution of the 
Hubble parameter. In fact, the likelihood posterior distributions of the parameters, and rj, are each 
separately compatible with vanishing conversion, even at la. The implication for this lack of reducibility 
is that the conversion can not be taken as a simple subleading effect, once it is included, the conversion 
must be non-vanishing to reach a best fit closer to the one obtained with the standard model. 

Despite the loss of direct compatibility among the standard parameters when comparing the two 



9 



500 



100 



0.01 



0.02 



0.05 



0.10 



0.20 



k{h Mpc"^) 



Figure 3: Matter power spectra for different models containing the best fit for WMAP+SDSS+SNIa. 
ACDM (magenta continuous line) and the model with dark matter conversion (black dashed). SDSS 
data are shown in blue points, but only diagonal elements of the covariance matrix. 



models, it is interesting to note how some of them can be relaxed or tightened, such as the dark 
energy equation of state parameter, the cold dark matter density or the neutrino masses. Although the 
parallel between best fits is the main criteria to qualify an alternative model, we show the relaxation of 
some parameters since it could be useful to explain some isolated effects related to specific parameters 
detached from the general model. The relaxation of neutrino masses bounds is interesting in the case 
one wants to add extra neutrino sterile families with masses at eV scale, such as suggested by anomalies 
in neutrino oscillations [39^ 00], although the allowed mass bound is still far from the scale of a few 
eV required by these oscillations anomalies. The conversion of cold into hot dark matter decreases the 
dark energy equation of state parameter to lower values, to compensate for the creation of matter whose 
density dilutes with a higher order of the scale factor. This effect is compatible with what is expected 
when the neutrino masses are increased, which also increases the amount of hot dark matter. 

The last and most interesting side effect of our conversion model would be the smaller value for the 
cold dark matter density at late times, specially in the central cores where it is expected to happen 
the conversion. The decrease in density is desirable given the incompatibility between the cusp galactic 
centers predicted by n-body simulations and the core profile observed in all surveys [41]. For instance, 
in the figure [5] we show the likelihood posterior distribution for all the standard model parameters, 
primitive and derived ones, with dark matter conversion, and the distribution for ACDM. It can be 
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Figure 4: Correlation of the parameters a, and ry at 68% CL and 95% CL. All others parameters 
are marginalized over their total mean values. The confidence regions correspond to the fits obtained 
with the following data sets: WMAP+SNIa (cyan), WMAP+SDSS (red) and WMAP+SDSS+SNIa 
(magenta) . 



seen that the largest deviation on primitive parameters happens for the total amount of dark matter 
density f^c- Given that the conversion decreases the amount of cold dark matter at late times, it is not 
unexpected that the fitting procedure found a best fit with small density today for cold dark matter 
while keeping its early time density compatible with the found in ACDM scenario. Such concentrated 
effect in cold dark matter density is not accidental, since the conversion model was built to diminish 
the dark matter density at small scales without affecting the nice agreement with large scale structures. 
The deviations from the other primitive parameters with respect to the scenario with no conversion are 
generally small and with hardly observable effects. 

6 Conclusions 

We have tested with large scale structure data the hypothesis of an environmental conversion of cold 
dark matter into hot dark matter that took place at late times in small scales. The effect was inspired 
in astrophysical processes that convert gravity binding energy in flow of relativistic particles. The main 
motivations to develop such phenomenological model are the particle models with composed dark matter 
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Figure 5: Posterior normalized distributions for all the standard model parameters for ACDM (dashed 
line) and for the model with dark matter conversion (solid line) . All others parameters are marginalized 
over their total mean values. 



candidates and the hints of dark radiation presence in different stages of the Universe evolution. One 
possible application is to address the discrepancy between the cold dark matter model and the very 
small scales observations, namely the core-cusp and the missing satellites issues. 

The results show that the pure cold dark matter is still a better option, although only marginally. The 
compatibility of the conversion with large scale structure was expected, since, aware that no conversion 
could take place in early times and in very large scales, we constrained the effect to happen at small 
scales and allowed to be triggered at late times, avoiding the most rigorous bounds in the data. In this 
way, the mechanism presented here works as a proposal for a model, already delimited by the available 
data at large scales. The main effects happen at very small scales and only observations at these scales 
will settle the complete portrait of the dark matter component. The results obtained in this work could 
hopefully help to direct the theoretical effort to address some issues at small scales when observational 
data will be available. The good agreement obtained with present data is an indication that this is not 
only valid, but a suitable solution. 

A possible extension of this work would be to perform the n-body simulation for some galaxy 
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clusters and put the conversion in the high density halos. Only accessing the non-linear regime of 
perturbation, one could explore the full capacity of such environmental mechanism of dark radiation 
production. Although the core-cusp and missing satellites issues are not the only observational hint for 
dark radiation, it is certainly the most interesting one for this model, given the factual relevance and 
for the possibility to be tested more directly with data. 
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